PROCEEDINGS 

OF SCIENCE 



in 
o 
o 

(N 

<D 

00 
(N 



> 
OS 

o 
in 
o 

D 



Real-time dynamics of proton decay 



Dmitri Grigoriev* 

Dept. of Mathematical Physics, NUI Maynooth, Co. Kildare, Ireland 
E-mail: |dima@thphys . may . ie 



Substituting Skyrmion for nucleon, one can potentially see — in real time — how the monopole 
is catalysing the proton (or neutron) decay, and even obtain a plausible estimate for catalysis 
cross-section. Here we discuss the key aspects of a practical implementation of such approach 
and demonstrate how one can overcome the main technical problems: Gauss constraint violation 
and reflections at the boundaries. 
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1. Introduction 

The monopole catalysis of proton decay [[[], is a complicated process involving several 
different physical scales. Although the physics of how fermions interact with the monopole is 
rather well understood, the net proton decay cross-section remains unknown, mainly because the 
proton decay actually happens via a combination of quite different physical mechanisms. At first, 
spatially separated monopole and proton interact at relatively large distances; then, the monopole 
starts to interact with individual quark states inside the proton, whose internal structure is defined 
by non-perturbative QCD phenomena. 

An interesting possibility to study both kinds of processes is provided [Q] by the Skyrme theory 
[Qj. However, even though the monopole catalysis of the Skyrmion decay is essentially a classical 
process, its cross-section is also not known. Recently it has been suggested ||] to study this 
process in a system of the Skyrmion interacting with the 't Hooft-Polyakov monopole [^, Hj. Since 
the latter is a spatially extended object, this makes the model free from singularities and greatly 
simplifies its numerical study. In the papers []|, ^], the spherically-symmetric monopole-Skyrmion 
geometry has been comprehensively studied in quasistatic limit. In particular, the decay path of 
the Skyrmion has been identified, and it has been shown that nothing prevents the Skyrmion from 
promptly losing almost all of its energy in the monopole background. 

The natural next step is to study the classical dynamics of the model suggested in H]. In 
this talk we discuss the main technical aspects of such a study, while the physical results will be 
published elsewhere []Tot]. In particular, we concentrate on a problem of gauge invariance violation 
in real-time classical simulations. 



2. The model 

The model 1 consists of the Skyrme field U coupled to two gauge fields (group SU (2)l) 
and Bp (group SU (2)r). The Higgs sector has two Georgi-Glashow Higgs fields <t>A and which 
break the SU (2)l and SU (2)r symmetries down to U(1)l and U(1)r, respectively, and a third Higgs 
field *F which breaks the axial subgroup of the remaining U(1)l x U(1)r. These requirements 
fix the scalar field group representations completely: U £ SU (2) , U, 1 ^ : (2,2), <I>a : (3,1), 
<S> B : (1,3). 

The action for the model is 
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For a detailed discussion of the structure and the properties of the model see the Refs. [pL n]. 
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Here F^ v and G^ v are the field strengths of A^ and B^, the fields <I>a and <J># are 2 x 2 matrices 
from the algebras of SU (2)l and SU (2)r, respectively, 
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(2.2) 



and F n and e are the pion decay constant and the Skyrme constant. The new Higgs field *F is a 
2x2 complex matrix with covariant derivative identical to that of the Skyrme field, 



Finally, the Wess-Zumino-Witten term Twzw can be safely ignored ^ in what follows. 

The most general spherically symmetric Ansatz consistent with the symmetries []5| ||] of the 



action (2.1) is 
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<3>A = $>b = ivh(r)x ■ o 
U = Icos/(r) + ix-5sin/(r) , »F = C,[ihi(r)x.o + h 2 (r)1] (2.3) 

where x is the unit radius-vector. 

Substituting the ansatz (2.3) into ( |Q| ) and using the Ao = Bo = gauge to get the proper 
Hamilton real-time evolution, one obtains spherically-symmetric action which remains invariant to 
the following static gauge transformations which are the remains of the axial t/(l): 

{h 2 + ih x )e ia ^ 
{a y +ia 2 )e ia[ - r) 
a?, + rd r a(r) 

f+a(r) (2.4) 

(note that h 2 + ih\ rotates identically to e'f , once the fields *F and U are in the same group represen- 
tation). This remaining gauge invariance becomes a major obstacle in real-time simulations (see 
below). 

The classical equations of motion can be obtained from the spherically-symmetric action by 
the standard variational procedure. Of particular importance for what follows is the Gauss con- 
straint corresponding to variation over ao- 
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The constraint (|2j) is preserved by the equations of motion as long as their gauge invariance is not 
violated. 

3. Numerical study 



In the Ao = Bq = gauge the classical equations of motion obtainable from the action ( |2.lD 
represent a system of 2 nd -order hyperbolic partial differential equations in (1+1) dimensions which 
can be solved numerically by e.g. a staggered leapfrog scheme. Choosing different initial config- 
urations, one can study several physically interesting cases of the Skyrmion decay. For example, 
starting from a normal vacuum Skyrme solution (so-called "bare" Skyrmion, see ^ ^]), one can 
simulate the decay of a physical nucleon. By choosing the initial profile for f(r) corresponding 
to the "thin" Skyrmion (a static solution [Q, ^] with unit baryon charge and minimal energy in the 
monopole background), it becomes possible to check [ |l0| ] how close is its actual decay path to the 
quasistatic paths found in the Refs.[]|, gj. However, any such study is facing two serious difficulties 
related to the lack of gauge-invariant lattice action for the model and the need for fully absorbing 
boundary conditions. 

3.1 Lattice action 

To preserve the physical structure of continuum theory, its symmetries must be retained in 
the discretised lattice version, and, in particular, the symmetries to local gauge transformations. 
Unfortunately, we are not aware of any gauge-invariant lattice action for gauged Skyrme mod- 
els, including the model used in the present study. The use of gauge-noninvariant lattice action, 
generally speaking, results in domination of nonphysical states in generated field configurations. 
However, in real-time classical evolution the mixing of physical and nonphysical states can be 
rather well controlled. 

Since the only reason for gauge invariance violation are the discretisation errors, for reason- 
ably small lattice spacing the magnitude of the violation remains rather small. Assuming that the 



evolution starts from a physical state which obeys the Gauss constraint ( |2.5[ ), the system initially 
stays close to the physical manifold and can be projected back onto it after certain amount of time. 
Physically, the projection procedure involves the minimisation of the quantity 

Q= fdxA 2 (3.1) 



where A is the constraint value (the left-hand side of (|2.5D). Generally speaking, any minimisation 
technique can be used here, with no specific approach favoured by physical reasons as long as the 
accumulated constraint violation remains small. In the present study we use a variation [Q] of the 
steepest descent method, also called sometimes as the Langevin cooling. 

It is worth noting that in real-time classical simulations the projection procedure is inevitable 
[^] even with a perfectly gauge-invariant lattice action. The reason is that the gauge invariance will 
inevitably be broken by computer roundoff errors. Even though the latter remain tiny (0(1O~ 15 ) 
for the standard double precision), the constraint violation accumulates in the course of evolution, 
which results in exponential divergency of constraints. To remain close to the physical manifold 
for indefinite period of time, one still has to periodically project the system onto it, although for 
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Figure 1: Constraint violation inevitably breaks stability of the numerical evolution, with the Gauss con- 
straint and total energy exploding simultaneously after the former exceeds a certain critical level. Of course, 
the numerical scheme also violates the energy conservation (as demonstrated by the small oscillations of en- 
ergy), but in non-gauge models this doesn't break the stability. Notations: K\ =F%/8v 2 , K2 =g 2 /64e 2 ; 
for dimensionless units used see Ref.n6|. 



gauge-invariant lattice schemes the constraint violation can be maintained at dramatically lower 
level. 

Failure to keep constraint violation under control inevitably leads to an explosion of the nu- 
merical scheme, whether a gauge-invariant one or not. For the model at hand this is demonstrated 
at the Figure 1, where the system energy is plotted against the time along with the magnitude of 
the volume-averaged constraint value < A > vo i- 

3.2 Absorption at the boundaries 

In physical case, the energy of the decaying proton is released to the spatial infinity. In nu- 
merical studies, this energy should be either absorbed at the boundaries or allowed to pass through 
the boundaries. For massless waves this can be easily implemented by the standard outgoing-wave 
boundary conditions. However, in more realistic cases, including the present study, it turns out to 
be very difficult to achieve full absorption at the boundaries, simply because any modification of 
the theory near the boundaries introduces spatial inhomogeneities which scatter the waves back 
to centre. Fortunately, in (l+l)-dimensional case it is still possible to overcome the problem by 
brute force, taking a sufficiently large physical volume, so that the reflected wave wouldn't come 
back to the origin before the end of the run, but this obviously cannot be done for higher spatial 
dimensions. 
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4. Conclusions 



Despite the problems discussed above, the monopole catalysis of Skyrmion decay can be suc- 
cessively studied [[K]] in the spherically-symmetric case where one can determine the decay time 
of the Skyrmion when it is exactly overlapped with the monopole. In two spatial dimensions, one 
can study Skyrmion-monopole collisions at zero impact parameter; however, the catalysis cross- 
section can be established only by studying full-3D geometry. This is a major challenge, taking 



into account that the Skyrme dynamics is known Q12[]-[|14|] to become unstable already in 2D simu 



lations, and a number of similar problems such as full numerical study of Skyrmion-anti-Skyrmion 
annihilation remain so far unsolved. 
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